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Abstract 

This report presents the results of tests at NASA Lewis 
to evaluate three icing scaling relationships or "laws" for 
an unheated model. The laws were LWC x time = con- 
stant, one proposed by a Swedish -Russian group and one 
used at ONERA in France. Icing tests were performed in 
the NASA Lewis Icing Research Tunnel (IRT) with 
cylinders ranging from 2.5- to 15.2-cm diameter. Refer- 
ence conditions were chosen to provide rime, mixed and 
glaze ice. Scaled conditions were tested for several sce- 
narios of size and velocity scaling, and the resulting ice 
shapes compared. For rime-ice conditions, all three of 
the scaling laws provided scaled ice shapes which closely 
matched reference ice shapes. For mixed ice and for 
glaze ice none of the scaling laws produced consistently 
good simulation of the reference ice shapes. Explana- 
tions for the observed results are proposed, and scaling 
issues requiring further study are identified. 

Nomenclature 

a Constant, dimensionless 

A Accumulation parameter, dimensionless 

C 

b Relative heat factor, dimensionless 

c Characteristic model length, cm 

c Specific heat, cal/g K 

h Convective film heat-transfer coefficient, 

c 2 
cal/secm K 

k Thermal conductivity, cal/sec m K 

LWC Liquid-water content, g/m 3 

M Mach number, dimensionless 

n Freezing fraction, dimensionless 

Nu Nusselt number, dimensionless 

2 

p Ambient static pressure, nt/m 

2 

p Ambient vapor pressure of water, nt/m 

Re Reynolds number, dimensionless 

V Airspeed, m/s 

t Ambient static temperature, °C 

/ Ambient total temperature, °C 

We Weber number, dimensionless 


/? o Droplet catch efficiency at stagnation point, 
dimensionless 

8 Droplet median volume diameter, pm 

p Density, g/cm 3 

p Viscosity, g/cm s 

a Surface tension of water against air, dyne/ cm 

w/a 

x Icing time, min 

Subscripts: 
a Air 

i Ice 

r Reference size and conditions 

s Scale size and conditions 

w Water 

Introduction 

This report presents the initial findings of some recent 
tests performed at NASA Lewis Research Center to 
investigate the validity of three icing scaling laws for an 
unheated model. 

In both wind tunnel and flight testing the researcher is 
often faced with facility or meteorological limitations 
which prevent testing at desired conditions. For tunnel 
tests, size restrictions rarely permit the testing of full-size 
models, and icing facilities can operate over only a 
limited range of air speeds and water-spray conditions. 
Flight testing eliminates the need to scale size and air- 
speed, but because cloud conditions cannot be controlled 
test-condition scaling is generally required. There is 
therefore a crucial need for reliable techniques to permit 
the scaling of size and test conditions in such a way that 
ice shapes resulting from scaled tests adequately simulate 
reference (full-size) results. 

A number of scaling laws have been derived by re- 
searchers over the past 30 years or so (see refs. 1-6.) 
Each of these laws consists of a set of relatively simple 
equations which relate scale test conditions to the refer- 
ence conditions when scale size and at least one other 
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scale variable (usually airspeed) have been defined. 

For this study, only three of the several published scaling 
laws were considered. They were: LWC x time = con- 
stant, the Swedish -Russian law 4 , and the French law 5 . 
Tests were performed in the Lewis Icing Research 
Tunnel (IRT) with cylinders as the test models. To eval- 
uate LWC x time = constant, several tests were conducted 
with different liquid-water contents and spray times but 
with cylinder size and all other test conditions the same 
from test to test. Liquid-water contents ranged from .4 to 
.7 g/m 3 for one series of tests and from .8 to 1 .3 g/m 3 for 
others. For the Swedish -Russian and French laws, the 
cylinder size and tunnel and cloud conditions were se- 
lected to define reference conditions and the scaled cylin- 
der size and airspeed chosen. Each law was then applied 
to determine the remaining scaled test conditions. 
Reference conditions included cylinder diameters of 1 5.2 
and 5. 1 cm (6 and 2 in), total temperatures of -26 to -8°C 
(- 1 5 to 1 8°F), velocities of 76 to 94 m/s ( 1 70 to 2 1 0 mph), 
median volume droplet diameters of 28 to 30 pm, liquid- 
water contents of .6 to 1 .3 g/m 3 , and spray times of 7.8 to 
18.3 min. These test conditions resulted in rime, mixed 
and glaze ice. Scaling scenarios tested included scale-to- 
reference size ratios of 1 : 1 and 1 :2 and scale-to-reference 
velocity ratios ranging from 1:1 to 1:2. The criteria for 
successful scaling were that the scaled ice shape should 
simulate the reference shape and that the amount of ice 
accreted in the scaled test, corrected for the model size, 
should match the reference quantity. 

In this repeal, the scaling laws will be described and ex- 
perimental ice shapes for reference and scaled conditions 
compared. Explanations for the observed results will be 
proposed, and scaling issues requiring further study will 
be identified. 

Considerations in the Development of Scaling Laws 

In an attempt to analyze the icing scaling problem rigor- 
ously, Bilanin 7 identified 18 dimensionless groups 
which affect the normalized ice thickness. A considera- 
tion of three of these similarity parameters, the Mach 
number, the Reynolds number and the Weber number, il- 
lustrates the fundamental difficulty of icing scaling anal- 


ysis: 

V 2 

M = (1) 
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Re ( 2 ) 
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V 2 cp 
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^w/a 

Because the temperatures of interest to icing extend over 
a very narrow range on the absolute temperature scale, 
T , c , u and a are effectively constant. Thus, to 

a jva *a w/a 

match the Mach number between scale and reference 
cases requires 

V - V (4) 

s r 

The altitudes at which icing occurs are generally low 
enough that test facility total pressures are close to flight 
icing-encounter total pressures. In this case, scale and 
reference p are nearly equal. Then to match the 
Reynolds number requires 

V = C ~V (5) 

s r 

C 

% 

and to match the Weber number requires 

.5 

y = y (6) 

s .5 r 
C 

s 

Clearly, it is impossible to satisfy similarity of Mach, 
Reynolds and Weber numbers simultaneously, except for 
the special case in which both scale size and velocity 
equal their respective reference values. The Mach num- 
ber can be neglected because at the relatively low speeds 
of interest to icing, compressibility effects are small. 
Even with this simplification, however, equations (5) 
and (6) are inconsistent, and each requires scale veloci- 
ties greater than the reference when the scale size is less 
than the reference. Because test velocities greater than 
typical flight airspeeds are not generally achievable in 
icing test facilities, equations (5) and (6) impose imprac- 
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tical restrictions. 

This dilemma is overcome in the existing analyses by ar- 
guing that not only are compressibility effects negligible, 
but viscous effects are as well, so that the scale and refer- 
ence Reynolds numbers need not be matched. The ratio- 
nale for this additional simplification has been that icing 
primarily occurs at the leading edge where the boundary 
layer is still thin 1 * * * * 6 . In addition, conventional scaling 
analyses make no attempt to match the Weber number; 
thus, they have assumed that water surface tension has 
no influence on ice growth. Instead of attempting to sat- 
isfy equations (1) - (3), existing scaling methods assume 
that the scale and reference cases will have similarity of 
geometry and of both pressure and velocity flowfields 
around the model. The scaling equations are then de- 
rived by applying the continuity (water-catch analysis), 
momentum (droplet-trajectory analysis) and energy (sur- 
face heat-balance analysis) equations to the ice-accretion 
process at the model stagnation line. The basic deriva- 
tions cm which most of the various laws are based has 
been given by Ruff 6 . 

Scaling Laws Tested 

The sets of equations used for each of the three laws con- 
sidered in this study will be given here. These equations 
have been incorporated into a NASA Lewis computer 
code to permit ready determination of scale conditions 
for different test cases. 

1 . Constant Product of Liquid- Water Content and Time 
This law considers the special case for which only scale 
liquid-water content and time are permitted to differ 
from the reference values. It is is the most restrictive of 
the three laws tested and is based on the least analysis. 

The analysis includes only water-catch considerations to 

insure that scale and reference ice accumulations are 
matched. However, this law also matches the droplet tra- 
jectory by requiring that the scale model size, airspeed, 
pressure and water droplet size be the same as the refer- 
ence values. The scale and reference temperature are re- 

quired to have the same value, but this is not sufficient to 

insure that the surface heat balances will match. Because 

of the restrictions in test conditions, this law also coinci- 


dentally satisfies similarity of Mach number, Reynolds 
number and Weber number; it is the only scaling law 
considered here to do so. 

The water-catch analysis identifies a term called the ac- 
cumulation parameter; it is defined as 

LWC V x 

A (7) 

c 

cp. 

The scale and reference accumulation parameters must 
be matched along with the catch efficiencies, to in- 
sure equal normalized ice thicknesses. For this law, the 
scale and reference values of model size, c, airspeed, V y 
and median volume droplet diameter, S y are matched; 
thus, the catch efficiency, which is a function of all three 
parameters, is also matched. Therefore, equation (7) 
shows that in order for the scale and reference normal- 
ized ice thicknesses to be the same, the product of liquid- 
water content and time must be the same. The complete 
set of equations to be satisfied is then 


c = c 

(8) 
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If the reference conditions, c , V , p , S , ( LWC ) , t and x 

r r r r r r r 

are known and a value for the scale liquid- water content, 
C LWC ) , is chosen, then all the scale parameters can be 
determined from eqs. (8-14). 

2 . Swedish -Russian This law was developed by a joint 
Swedish-Soviet working group on flight safety 4 . In ad- 
dition to the water-catch considerations of the previous 
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law, it includes a droplet-trajectory analysis. Energy- 
balance similarity is approximately satisfied by setting 
the scale temperature and liquid-water content equal to 
their respective reference values. The user selects the 
scale size and airspeed. The static pressure in most facil- 
ities, including the IRT, cannot be controlled, so it is de- 
termined by airspeed and ambient pressure. The 
complete scale conditions are then given in the following 
equations: 


= [selected by user] 

V = [selected by user] 
p ** [fixed by conditions] 
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3. French The French scaling law was developed by re- 
searchers at ONERA 5 . This law includes an analysis for 
surface energy balance in addition to the analyses for 
water catch and droplet trajectory. The energy-balance 

g 

analysis is based on the work of Messinger . Unlike the 
other laws considered here, no scale parameters are re- 
quired to equal their reference values. As used in this 
study, the French scaling law is: 
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Because the ambient vapor pressure of water, p , is a 
function of temperature, equation (28) must be solved it- 
eratively to find the air temperature of the scale case. 


The development of equation (28) from the Messinger 
analysis requires that the scale values of freezing frac- 
tion, n, and relative heat factor, b, be equal to their re- 
spective reference values. The freezing fraction is the 
fraction of water arriving at a given location on the 
model which freezes at that location. The relative heat 
factor was defined by Tribus 9 as 


b = 


LWC V Be 

r r 0 p,w 


h 

c 


(29) 


Here, c is the specific heat of water at the surface 

p,w 

temperature of the model, and h is the surface 
convective heat transfer coefficient, h can be found 

c 

from the Nusselt number, Nu, using the usual definition. 


h =Nuk fc (30) 

c a 

Equation (26) results from equating the scale and refer- 
ence relative heat factors. It can be seen from the expo- 
nents in this equation that the Nussselt number and 
Reynolds number. Re, are assumed to be related as 

Nu=aRe ' 8 (31) 
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In the studies reported here, a value of .033 was used for 
the constant a. In equation (3 1 ) the Reynolds number is 
based on the airspeed, V y and the characteristic 
dimension of the model, c: 

Re^Vcpffi (32) 

a a 

The scale static temperature found from equation (28) is 
insensitive to the value of b from equation (29). Thus, al- 
though the droplet catch efficiency, that appears in 
equation (29) can be determined by the method given by 
Langmuir and Blodgett 10 , for this study it was suffi- 
ciently accurate to use a value of fi = .5. 

The derivation of equation (28) assumes that the surface 
is at the freezing temperature of water. This assumption 
is only valid if the freezing fraction, n, is less than one. 
Because of this assumption about the surface tempera- 
ture, equation (28) is not strictly applicable to rime-ice 
conditions. However, for rime conditions water freezes 
immediately on impact, and the ice shape is insensitive 
to the ambient temperature. Thus, the temperature found 
from equation (28) is adequate for scaling as long as the 
scaled ice is also rime. In this study, equation (28) was 
used for rime as well as for mixed and glaze ice. 

Description of Experiment 

NASA Lewis Icing Research TunneL The experiments 
were performed in the NASA Lewis Icing Research 
Tunnel (IRT) shown in figure l 11 . The IRT has a test 
section width of 2.74 m (9 ft) and a height of 1.83 m (6 
ft.) With an empty test section, it is capable of operation 
at test-section velocities up to 130 m/s (300 mph). A 
refrigeration system permits accurate control of the test- 
section temperature from -40 to 0°C (-40 to 32°F). A 
water- spray system 12 with 8 spray bars containing ap- 
proximately 100 spray nozzles provides the ability to 
control test-section liquid-water content from .2 to 
3 g/m 3 and droplet median volume diameters from 15 to 
40 Jim. 

Two sets of spray nozzles, known as the mod-1 and stan- 
dard nozzles, are used in the IRT to provide different 
ranges of liquid-water content and droplet size . Both 


nozzles were used in these tests. Some spray conditions 
can be obtained with either nozzle set, and studies have 
shown that the ice shape for a given set of conditions is 
independent of the nozzle type 1 3 . 

Scaling Test Hardware. Ice accretion was measured on 
hollow circular aluminum cylinders with 1 5.2-, 7.6-, 5. 1 - 
and 2.5-cm (6-, 3-, 2- and 1-in) diameters. Each cylinder 
was mounted vertically in the center of the test section as 
shown in figure 2. As seen in figure 2(a), the 15.2-, 7.6- 
and 5. 1 -cm (6-, 3- and 2-in) cylinders were installed with 
a short mounting cuff at each end. The cuffs were at- 
tached to aluminum angles mounted on the ceiling and 
floor of the tunnel test section. These cylinders extended 
over the full height of the test section. The 2. 5-can ( 1 -in) 
cylinder was installed as shown in Figure 2(b). It hung 
from a sleeve mounted to a plate attached to the test-sec- 
tion ceiling, and it extended just below the center of the 
tunnel. All cylinders were marked with a line around 
their circumference to identify the middle of the test sec- 
tion. All ice tracings were made at this location. 

T est Procedure. Tests were performed by first 
establishing the tunnel velocity and total temperature. 
The pressures of the spray-bar water and atomizing air, 
which determine LWC and 8, were then set, and when 
tunnel conditions had stabilized, the water spray was 
initiated. To account for the start-up transient, spray 
times for the tests were increased by 1 2 s from the values 
determined from the scaling laws. This addition was 
never more than 4 percent of the total spray time, and it 
will not be included in the times reported here. When the 
prescribed spray period was completed, the spray was 
shut off and the tunnel brought to idle to permit person- 
nel entry into the test section. After the ice shape had 
been recorded, the ice was cleaned from the model and 
the procedure repeated for the next spray condition 

The ice shape was recorded manually for each test. A 
heated aluminum block with a semicircular cut-out 
matching the cylinder was used to melt a thin cut into the 
ice normal to the cylinder axis at the center of the test 
section. A cardboard template, also with a semicircular 
cutout to match the cylinder, was placed in this gap, and 
the ice shape was traced onto the cardboard template. 
The tracing was later digitized for computer storage. 
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Results 

Test Repeatability Evaluations in this study were based 
on judging differences in ice shape from one run to an- 
other. Thus it was also necessary to establish, for a 
benchmark, the typical variation in ice shapes from run 
to run when identical tunnel conditions were specified. 
Ice accretes on tunnel components during the course of 
each day’s testing, possibly altering the tunnel flow and 
spray characteristics. It might also be possible for the 
tunnel or spray bar characteristics to change for a variety 
of reasons over longer periods of time. To address these 
possibilities, extensive repeatability tests have been 
performed in the IRT by Shin and Bond 13 . Their tests 
showed that repeatability of ice shapes cm airfoils was 
generally excellent. In view of these results, it was not 
necessary to perform additional comprehensive repeabil- 
ity testing; however, to achieve the purposes of the pre- 
sent study, several sets of test conditions were repeated 
from time to time. 

Figure 3 presents typical data for repeatability tests 
performed on the 1 5.2-cm- (6-in-) diameter cylinder (fig. 
3(a)) and the 5.1 -cm (2-in) cylinder (fig. 3(b).) The 
1 5.2-cm (6-in) cylinder tests reported in figure 3(a) show 
that the amount of ice accreted and the lower horn size 
were consistent from run to run . The upper horn size and 
location varied somewhat, however. 

Figure 3(b) shows results for the same operating 
conditions as 3(a) but for a 5.1 -cm- (2-in-) diameter 
cylinder. Only two tests were made for this cylinder at 
these conditions, and very little variation in the ice 
shapes was observed. Additional repeatability tests were 
also made for other cylinders at various conditions; the 
results shown here are representative of the ice-shape 
variation observed. 

Scaling With Constant L IT^xTime Figure 4 gives three 
sets of results for the LWC x time = constant law. All 
tests were made with a tunnel airspeed of 94 m/s (210 
mph). 

Figure 4(a) compares results on a 1 5.2-cm (6-in) cylinder 
at a total temperature of -26°C (-14°F) and a droplet me- 
dian volume diameter of 30 pm. Three liquid-water 


contents were tested, .8, 1.0 and 1.3 g/m 3 , and therespec- 
tive spray times were 12.7, 10. 1 and 7.8 min. For these 
rime-ice conditions, both the quantity of ice accreted and 
the ice shape were preserved for all three tests, showing 
that this scaling law applies for these conditions. 

Ice shapes on the 5.1 -cm- (2-in-) diameter cylinder for a 
somewhat higher temperature and lower liquid-water 
contents are shown in figure 4(b). At the lowest liquid- 
water content shown, .4 g/m 3 , the ice shape suggested 
that little water runback occurred. When the liquid- wa- 
ter content was increased to .5 g/m , the higher water 
loading reduced the freezing fraction, and water runback 
from the stagnation region increased. As a result, a 
small valley formed in the ice at the stagnation line with 
stubby horns adjacent to it. Finally, a liquid-water con- 
tent of .7 g/m 3 produced a still deeper valley and even 
more pronounced horns. Although total ice accumula- 
tion appeared to be constant for these three tests, the ice 
shape was not. 

Figure 4(c) shows ice shapes on the 5.1 -cm (2-in) 
cylinder at the same temperature as figure 4(b) but with a 
different drop size and with liquid-water contents of .8, 1 
and 1.3 g/m 3 . These conditions resulted in fully glaze ice 
shapes. The ice thickness near the stagnation line was 
similar for all three tests, and the total ice accumulation 
appeared to be the same. However, the angle between the 
horns increased significantly when liquid-water content 
increased from .8 to 1 g/m and from 1 to 1.3 g/m . As 
with the mixed-ice conditions of figure 4(b), these results 
illustrate how liquid- water content can have a significant 
effect on ice shape. 

LWC x time = constant satisfies Mach-, Reynolds- and 
Weber-number similarity, and both water catch and 
droplet trajectory are matched. Although it's not clear if 
all these features are required, application of this law 
produced a constant amount of ice for a range of liquid- 
water contents for rime, mixed and glaze ice. However, 
except for fully rime conditions, ice shapes could not be 
scaled with this law. The reason is that the liquid-water 
content can have a significant effect on the surface heat 
balance, and this law does not match surface heat bal- 
ance terms between scale and reference conditions. 
These results show that it is not sufficient to scale icing 
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encounters by satisfying similarity in Mach number, 
Reynolds number, Weber number, droplet trajectory and 
water catch alone. It is also necessary that the surface 
heat balance for the scale test be matched to the refer- 
ence. 

Scaling Using Swedish-Russian and French Laws 

These two laws do not satisfy similarity in Mach number, 
Reynolds number and Weber number, but each includes 
some consideration of the leading-edge heat balance. 
Results of tests to evaluate the Swedish-Russian and 
French laws are shown in figures 5, 6 and 7. In each fig- 
ure, the ice shape for the scaled tests are compared to the 
reference shape. 

Figure 5 gives rime-ice results. In figure 5(a) the model 
size was scaled from 1 5.2 to 7 .6 cm (6 to 3 in) and the air- 
speed was scaled from 76 to 61 m/s (170 to 136 mph). 
The quantity of ice accreted and the shape of the scaled 
ice agreed very closely with the reference ice for both 
laws. 

Scaling the size from 5. 1 to 2.5 cm (2 to 1 in) and the air- 
speed from 94 to 47 m/s (210 to 105 mph) produced the 
rime ice shapes shown in figure 5(b). For this situation, 
the scaled ice shapes were a reasonable, although not 
perfect, match of the reference shape. It appears that the 
reference conditions resulted in slightly less freezing oc- 
curring along the stagnation line than occurred for either 
of the scaled conditions. However, it can be concluded 
that for rime conditions, both of these scaling laws pro- 
vided adequate scaling guidance. 

Figure 6 gives results at mixed-ice conditions. For both 
figures 6(a) and 6(b), the size was scaled from 5. 1 to 2.5 
cm (2 to 1 in). In figure 6(a), the scale and reference air- 
speeds were matched at 94 m/s (210 mph) while for fig- 
ure 6(b) airspeed was scaled from 76 to 61 m/s (170 to 
1 36 mph). For the scaling situations of both figures, the 
quantity of ice accreted seemed to match the reference 
quantity fairly well when the French law was used, while 
the Swedish-Russian law may have produced somewhat 
too much ice. 

For the cases presented in figures 6(a) and 6(b), the 
scaled ice shapes for both laws only approximately simu- 


lated the reference shapes. The reference shapes had 
deeper troughs and more pronounced horns than either 
of the scaled shapes, suggesting that the reference case 
had a lower freezing fraction at the leading edge. The 
Swedish-Russian law makes no attempt to match scale 
and reference freezing fraction; it satisfies the heat bal- 
ance by simply matching scale and reference static tem- 
perature and liquid-water content. As noted earlier, the 
French law does match the freezing fraction. However, 
it's not sufficient just to match the freezing fraction; the 
physical model on which it is based must also be correct. 

The freezing fraction depends in part on the value of the 
heat-transfer film coefficient. As noted earlier, the 
French law assumes the Nusselt number to be propor- 
tional to the Reynolds number raised to the .8 power. Ex- 
perimental airfoil data of Gelder and Lewis 14 showed 
that for turbulent boundary-layer flow, Nu oc Re’*, while 
for laminar boundary-layer flow (near the leading edge), 
Nu x Re' 5 . Gelder and Lewis based their Reynolds num- 
bers on the local velocity and their Nusselt and Reynolds 
numbers on the distance along the surface from the lead- 
ing edge. Van Fossen, et. al . 15 measured heat transfer on 
cylinders with and without artificial ice shapes for 
Reynolds numbers up to about 1 . 9x 1 0 5 . These Reynolds 
numbers were based on the cylinder diameter and free- 
stream velocity. At the leading edge the Nusselt number, 
also using cylinder diameter, varied with the Reynolds 
number raised to a power near .5 for low turbulence and 
unroughened model surfaces. In the tests reported here, 
the Reynolds number based on cylinder diameter and 
free-stream velocity (eq. (2)) varied from approximately 
10 5 to 10 6 Thus, for these conditions, the .5 power may 
be more appropriate than the .8 power used by the French 
scaling law. In its present form, the French scaling law 
may be more suitable for high-Reynolds-number appli- 
cations. It should be noted that the AEDC-law 6 analysis 
was based on a power of .5; thus, it may be a better choice 
for scaling at low Reynolds numbers. 

The effect of changing the Reynolds-number power from 
.8 to .5 in the French analysis can be seen by equating the 
scale and reference relative heat factors using Nu x Re' . 
Equation (26) then becomes 
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If equation (33) had been used instead of equation (26), 
LWC for the case of figure 6(a) would have been .85 in- 

stead of .7 g/m . For figure 6(b), LWC would have been 

3 * 

1.25 instead of .95 g/m . Higher liquid-water contents 
decrease the freezing fraction, and the lower freezing 
fractions would tend to reduce the leading-edge ice thick- 
ness and increase the horn size. This effect was illus- 
trated in figure 4(b) which showed ice shapes for mixed- 
ice conditions with three liquid-water contents. Thus, 
the change in Reynolds-n umber power from .8 to .5 is in 
the right direction to improve the agreement between the 
scaled and reference ice shapes, but additional test data 
are needed to confirm the validity of this change. High- 
Reynolds-number data to verify the existing French law 
and low-Reynolds-number data for the AEDC law are 
needed as well. 


Typical glaze-ice test results are presented in figure 7. 
Reference conditions for the tests of both figure 7(a) and 
(b) were identical except for the reference model size. 
Figure 7(a) presents results for scaling velocity from 94 
to 58 m/s (210 to 1 30 mph) with the cylinder diameter of 
15.2-cm (6-in) matched. For this situation, application 
of the scaling laws produced somewhat more ice than the 
reference test while the scale horns were less distinct 
than the reference. The reference conditions for figure 
7(a) were the subject of the repeatability study reported in 
figure 3(a). It is evident that the differences between the 
scaled and reference ice accretions in figure 7(a) are 
greater than the typical run-to-run variations shown in 
figure 3(a). 

Glaze results for the 5. 1-cm (2-in) cylinder are shown in 
figure 7(b). Both size and velocity were scaled in this 
case. Again, both scaling laws led to a significantly 
larger ice accretion than the reference conditions, and 
the prominent horns of the reference ice shape have been 
totally lost. The comments about the poorly-scaled 
shapes of figures 6(a) and (b) also apply here and need 
not be repeated. Possible reasons for the poor simulation 
of the total amount of ice accreted will be discussed fur- 
ther. 


Both the Swedish-Russian and French laws match the 

accumulation parameter, A , and the droplet catch effi- 

c 

ciency, ft 9 between scale and reference conditions. The 
close match of reference and scaled ice accretions for 
rime and mixed conditions verify that these effects have 
been properly addressed. Therefore, the excess quantity 
of scaled ice accreted for the glaze cases must have re- 
sulted from neglected phenomena which are of impor- 
tance only for glaze ice. 

One such phenomenon is droplet splashing. If droplet 
splashing has a significant effect on ice accretion for 
glaze ice, the Weber number would need to be matched 
between scale and reference tests . For the present tests 
the Weber number for the scale tests was lower than the 
reference, and the tendency for splashing in the scale 
tests was therefore less than in the reference. 

Another mechanism potentially leading to reduced ice 
accretion is shedding of liquid water from the sur- 
face 16,17 . The scaled Reynolds numbers for both figures 
7(a) and (b) were smaller than their respective reference 
values. Consequently, the surface water in the scaled 
tests would have experienced lower shear than if Re had 
been matched. 

The ability of the LWC x time = constant law, which sat- 
isfies similarity in both Weber and Reynolds numbers, to 
maintain constant quantities of ice provides additional 
evidence that these parameters may be too important to 
neglect for glaze-ice conditons. Studies of the effects of 
both the Weber and Reynolds numbers are needed to im- 
prove our understanding of the physics of ice accretion. 
A better physical understanding is essential not only for 
the development of more effective scaling methodologies 
but also for the evolution of more accurate ice-accretion 
codes. 


Concluding Remarks 

Three scaling laws, LWC x time = constant, the Swedish- 
Russian law and the French law, were used to scale from 
reference conditions representing rime, glaze and 
mixed-ice conditions for an unheated model. Tests of the 
Swedish-Russian and French laws included scaling of 
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size alone, velocity alone, and of both size and velocity. 

The three laws represent a progression in analytical thor- 
oughness with the LWC x time = constant law involving 
the amplest analysis. This law satisfies similarity of 
Mach number, Reynolds number and Weber number. It 
always provided accurate simulation of the amount of ice 
accreted. However, it includes no consideration of heat 
balance; consequently, for mixed or glaze-ice conditions 
it did not scale ice shapes adequately for the liquid-water 
contents tested (.4 to 1 . 3 g/m 3 ). This law can be expected 
to be valid in general only for rime conditions, for which 
water freezes instantly on impact; i.e., for situations for 
which the surface heat balance has no effect on the ice 
shape. 

The Swedish-Russian and French laws gave nearly the 
same results. For rime ice, both laws produced scaled ice 
shapes that closely matched the reference. For mixed 
and glaze ice, the scaled ice shapes were only approxi- 
mate simulations of the reference shapes, having less- 
pronounced horns. Because the form of the Nusselt 

Q 

number, Nu x Re , used in the French-law analysis is 
applicable to higher Reynolds numbers than were tested 
in this study, it was speculated that better results might 
be achieved with a modified form of this law using a 
Reynolds-number power of .5 instead. Additional tests 
at the Reynolds numbers tested here, 10 5 to 10 6 , are 
needed to evaluate the effect of the Reynolds-number 
power on scaling results, and tests at higher Reynolds 
numbers are needed to test the validity of the French scal- 
ing law in its present form. 

For glaze conditions not only was the shape incorrectly 
scaled, but the quantity of scaled ice accreted was too 
great when either the Swedish-Russian or French laws 
was used. Neither of these laws satisfies similarity of 
Reynolds or Weber number. Because the LWC x time = 
constant law satisfies similarity of Reynolds and Weber - 
number and also produced the correct amount of glaze 
ice, it was reasoned that features of the ice-accretion pro- 
cess which depend on the Reynolds and Weber numbers 
may need to be included in analyses to adequately scale 
for glaze-ice conditions. In particular, phenomena such 
as droplet splashing and liquid-water shedding were con- 
sidered by Bilanin 7,16 and Olsen 17 to be relevant to ice 


accretion. Tests are needed to establish their importance. 
It is not clear how a practical scaling methodology which 
includes Reynolds and Weber number matching could be 
developed; however, studies to increase our understand- 
ing of the ice-accretion process should lead not only to 
improved scaling methodologies but also to more accu- 
rate ice-accretion codes. For many situations it may not 
be possible to derive exact scaling procedures, but ap- 
proximations based on a good physical understanding 
and supported by accurate analytical tools should prove 
to be adequate to satisfy the needs of icing experiments. 
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Vent Door 



p M 130-m/s (300-mph) Test Section 

Room [— [mu] | l g3 x 2 ?4 m (6 x 9 ft) 

Figure 1. NASA Lewis Icing Research Tunnel. 



Test Cylinder Control Room Mounting Sleeve Test Cylinder 

(a) 5.1- to 15.2-cm (2- to 6-in) Cylinders. (b) 2.5-cm (1-in) Cylinder. 

Figure 2. Test Cylinder Installed in IRT Test Section. 





8-11-92 Run 6 

1-6-93 Run 2 

1-6-93 Run 12 


1-7-93 Run 13 

1-8-93 Run 14 


(a) Cylinder Diameter, 15.2 cm (6 in). (b) Cylinder Diameter, 5.1 cm (2 in). 

Figure 3. Typical Ice Shape Repeatability. Total Temperature, -8°C(18°F); Velocity, 94 m/s (210 mph); 

Droplet Median Volume Diameter, 30 pm; Liquid-Water Content, 1.3 g/m 3 ; Icing Spray Time, 7.8 min. 
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(a) Rime Ice; Cylinder Diam., 15.2 cm (6 in); (b) Mixed Ice; Cylinder Diam., 5.1 cm (2 in); Total Temp., 

Total Temp., -26°C (-14°F); Droplet Median -8°C (18°F); Droplet Median Volume Diam., 20 pm. 

Volume Diam., 30 pm. 
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(c) Glaze Ice; Cylinder Diam., 5. 1 cm (2 in); Total Temp., 

-8°C (18°F); Droplet Median Volume Diameter, 30 pm. 

Figure 4. Scaling With LWC x Time = Constant. Velocity, 94 m/s (210 mph); LWC x r, 10.15 g min/m 3 . 
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(a) Reference Cylinder Diain., 15.2 cm (6 in). (b) Reference Cylinder Diam., 5.1 cm (2 in). 

Figure 5. Scaling both Size and Velocity With Swedish-Russian and French Laws. Rime Conditions. 
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(a) Scaling Size With Velocity Matched. (b) Scaling Both Size and Velocity. 

Figure 6. soling With Swedish-Russian and French Laws. Mixed-Ice Conditions; Reference Cylinder Diameter, 5. 1 cm (2 in). 





c V 8 LWC t 

cm °C m/s |im g/m 3 min 

Reference 15.2 -8 94 30 1.3 7.8 

Swedish-Russian 15.2 -11 58 38 1.3 12.6 

French 15.2 -10 58 36 1.4 11.5 



Reference 5.1 -8 94 30 1.3 7.8 

Swedish-Russian 2.5 -11 58 27 1.3 6.3 

French 2.5 -10 58 23 1.6 5.0 

(b) Scaling Both Size and Velocity. 

Reference Cylinder Diameter, 5.1 cm (2 in). 


(a) Scaling Velocity With Size Matched. 

Reference Cylinder Diameter, 15.2 cm (6 in). 


Figure 7. Scaling With Swedish-Russian and French Laws. Glaze-Ice Conditions. 
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